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> '. Abstract 

lO \ Geological fault systems, as the San Andreas fault (SAF) in USA, constitute typical examples of 

'nI" ■ self-organizing systems in nature. In this paper, we have considered some geophysical properties 



of the SAF system to test the viability of the nonextensive models for earthquakes developed in 
[Phys. Rev. E 73, 026102, 2006]. To this end, we have used 6188 earthquakes events ranging in 
the magnitude interval 2 < m < 8 that were taken from the Network Earthquake International 

rH ' Center catalogs (NEIC, 2004-2006) and the Bulletin of the International Seismological Centre (ISC, 

O '■ 

C_> , 1964-2003). For values of the Tsallis nonextensive parameter q ~ 1.68, it is shown that the energy 

distribution function deduced in above reference provides an excellent fit to the NEIC and ISC 

'i_l ■ SAF data. 
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I. INTRODUCTION 

Earthquakes are among the most relevant paradigms of the so-called self-organized criti- 
cality, introduced by Bak, Tang and Wiesenfeld in Ref. (l|. In the context of fault systems, 
this represents a complex spatio-temporal phenomenon which are related with the deforma- 
tion and sudden rupture of some parts of the earth's crust driven by convective motion in the 
mantle, i.e., the radiation of energy in the form of seismic waves. In particular, fault system 
as the San Andreas fault in California (USA) is one of the first example of self-organizing 

n 

system in nature 2]. Despite the complexity of earthquakes, some of the known empirical 
laws are considerably simple as, for instance, the Omori law [3| for temporal distribution 



a 



for relationship between frequency and 



of aftershocks and the Gutenberg-Richter law 
magnitude. 

The physics of earthquakes has also been studied in the context of nonextensive for- 
.a>..s m, wKe. the ve,v «.. ,„ve.tiga«o. wa. doe by Abe . Re.. 030. Reeen.y, 
a very interesting model for earthquake dynamics related to Tsallis nonentensive framework 
has also been proposed in Ref. [lO|. More recently, such a model was revisited by consid- 
ering some new ingredients, namely, the correct definition for mean values in the context 
of Tsallis nonextensive statistics, and a new scale between the earthquake energy and the 
size of fragments, e oc r^; see Ref. |ll] for details. In reality, by using the standard method 
of entropy maximization a new energy distribution function (EDF) was deduced, differing 
considerably from the one obtained in Ref. [lO|. 







In this paper, we test the viability of our approach [Ijj by using data taken from the 
Network Earthquake International Center catalogs (NEIC, 2004-2006) and the Bulletin of 
the International Seismological Centre (ISC, 1964-2003) for one of the best studied major 
fault zone in the world, i.e., the San Andreas Fault (SAF). As is well known, the SAF is one 
of the world longest and most active geological faults, reaching ~15 kilometers deep with an 
extension of ~1200 km and being about 20 million years old. This fault forms the boundary 
between the North American and Pacific plates and is classified as a right lateral strike-slip 
fault (transform boundary), although its movement also involves comparable amounts of 
reverse slip jl2|- In the geological and geophysical context, a considerable number studies 
have been performed in order to better understand this complex system (see, e.g. 
references therein). 



m and 



II. NONEXTENSIVE DESCRIPTION FOR EARTHQUAKES: A MODEL 

In this Section, we recall the theoretical basis of our nonextensive earthquake model 



discussed in Ref. 



Q 



As well known, the Tsallis' framework generalizes the Botzmann- 
Gibbs statistics in what concerns the concept of entropy. Such formalism is based on the 
entropy formula given by 

S, = - j p\a)\Ti,p{a)da. (1) 

Here and hereafter, the Boltzmann constant is set equal to unity for the sake of simplicity, 
q is the nonextensive parameter and the g-logarithmic function above is defined by 

ln,p= ^P'''-^\ (p>0), (2) 

1 — g 

which recovers the standard Boltzmann-Gibbs entropy 

5*1 = — / phipd^p 



in the limit g ^ 1. In our model, p(cr) stands for the probability of finding a fragment 

of relative surface a, which is defined as a characteristic surface of the system. It is worth 

mentioning that most of the experimental evidence supporting Tsallis proposal are related to 

he power-law distribution associated with 5*^ descripition of the classical A^-body problem 

^'- nn 

The fundamental idea of the model developed in Refs. [10, |ll| consists in the fact that 



the space between faults is filled with the residues of the breakage of the tectonic plates. In 
this regard, it was studied the influence of the size distribution of fragments on the energy 
distribution of earthquakes. The theoretical motivation follows from the fragmentation 
phenomena in the context of the geophysics systems [la]- In this latter work, Englman et 
al. showed that the standard Botzmann-Gibbs formalism, although useful, cannot account 
for an important feature of fragmentation process, i.e., the presence of scaling in the size 
distribution of fragments, which is one of the main ingredients of our approach. Thus, a 
nonextensive formalism is not only justified in these earthquake models but also necessary 



^ For a complete and updated list of references on Tsallis' entropy see http.tsallis.cat.cbpf.br/biblio. htm. 



since the process of violent fractioning is very probably a nonextensive phenomenon, leading 
to long-range interactions among the parts of the object being fragmented (see, e.g., |l7[). 



In Ref. [Ul, however, differently from [10'], which assumes e ~ r, we used a new energy 
scale e ~ r^. Thus, the proportionality between the released relative energy e and r^ (r is 
the size of fragments) is now given by a — cTg = (e/a)^''^, where a scales with r^ and a (the 
proportionality constant between e and r^) has dimension of volumetric energy density. In 
particular, this new scale is in full agreement with the standard theory of seismic rupture, 
the well-known seismic moment scaling with rupture length (see, for instance, |l8|). 

Following the standard method of maximization of the Tsallis entropy under the con- 
strains of the normalization of the function and g-expectation value, we obtain (see jll| for 
details) 



:^-i^'^-- 



1-9 

(3) 



p{a) 

which corresponds to the area distribution for the fragments of the fault, and the EDF of 
earthquakes reads as follows 



p{e)de = —, (4) 

which has also a power-law form with C and C given by 

C = ^ and C- = -p<i^. (5) 

In the above expression, the energy probability is written as p{e) = n{e)/N, where n{e) 
corresponds to the number of earthquakes with energy e and N total number earthquakes. 



III. TESTING THE MODEL WITH EXPERIMENTAL DATA 

In order to test the viability of the new EDF above derived [Eq. (j^] we introduce the 
cumulative number of earthquakes, given by the integral |li| 



Ne> 



p{e)de, (6) 



or, equivalently, 

(7) 
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FIG. 1: The relative cumulative number of earthquakes [Eq. ©] as a function of the magnitude 
m. The data points correspond to 6188 earthquake events of San Andreas fault (lying in the 
interval 2 < ?7i < 8) from Network Earthquake International Center catalogs (NEIC, 2004-2006) 
and Bulletin of the International Seismological Centre (ISC, 1964-2003). The best-fit values for 
the parameters q and a are shown in the Panel. A summary of this and other recent analyses is 
also shown in Table I. 

where N^y is the number of earthquakes with energy larger than e. By considering m = 
3^^ loge (m stands for magnitude) it is possible to rewrite the above Eqs. as (g 7^ 1) 

'2-q 



log(Ar>J = \ogN + 



X log 



1- 



1-q 
2-q 



(8) 
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which, similarly to the modified Gutenberg- Ricther law (see, e.g., Refs. J2l| for more details), 
describes appropriately the energy distribution in a wider detectable range of magnitudes. 
Figure 1 shows the main results of the present analysis. The relative cumulative number 
of earthquakes (defined as Gm> = N^^/N) is shown as a function of the magnitude m. 



The data points, corresponding to earthquakes events (produced by the San Andreas fault) 
ranging in the interval 2 < m < 8 with 6188 events, were taken from the Network Earthquake 
International Center catalogs for preliminary determinations (NEIC, 2004-2006) and Bulletin 
of the International Seismological Centre (ISC, 1964-2003). Note that we have used only the 
data from strike-slip movement of the San Andreas fault. The best-fit parameters for this 
analysis are given by g = 1.68 and a = 4.4 x 10^°, which is very close to the values found 
for other geological faults at different parts of the globe. Note also that, similarly to our 
previous analyses, the nonextensive model of Ref. [ll|, represented by Eqs. Q-®, provide 
an excellent fit to the experimental data of the San Andreas fault. 

At this point we should compare the predictions of the model discussed above with the 
ones of Ref. [li|. We note that a similar analysis to the latter model would provide the 
very same value for the nonextensive parameter, i.e., q = 1.68, although the energy density 
a would differ by several orders of magnitude, a = 4.3 x 10^^. Naturally, this huge difference 
between the predictions of these nonextensive earthquake models could be used to check 
the viability of each model by comparing their predictions with experimental estimates for 
the energy density a. Note also that from most of the geophysical analyses performed so 
far values of g ~ 1.6 — 1.7 seem to be universal, in the sense that different data sets from 
different regions of the globe indicate a value for the nonextensive parameter lying in this 
interval. It is possible, however, that the predictions of the nonextensive earthquake models 
discussed in this paper be dependent on the fault mechanisms generating the earthquakes. 
A detailed analysis on the influence of the fault mechanisms (reverse, normal or strike-slip) 
on the predictions of these nonextensive models or, more specifically, on the values of the 
nonextensive parameter q will appear in a forthcoming communication 

IV. CONCLUSION 

We have tested the viability of the earthquake model developed in Ref. |11| from the best 
studied major fault zone in the world, i.e., the San Andreas fault. By using 6188 earthquake 
events (in the interval 2 < -m < 8) taken from the Network Earthquake International Center 
catalogs and Bulletin of the International Seismological Centre we have shown, in agreement 
with other similar analyses, that for values of the nonextensive parameter of the order of 



q = 1.6 — 1.7 and a ~ 10 , the model of Ref. |ll| provides an excellent fit to the strike-slip 



TABLE I: Limits to q and a 



Fault 


Ref. 


q 


a 


California - USA 


[10] 


1.65 


5.73 X 10-6 


Iberian Penisula - Spain 


[10] 


1.64 


3.37 X 10-6 


Andalucia - Spain 


M 


1.60 


3.0 X 10-6 


Samambaia - Brazil 


lUJ 


1.66 


1.8 X 10^° 


New Madrid - USA 


lUJ 


1.63 


1.2 X 10^° 


Anatolian - Turkey 


11^ 


1.71 


2.8 X 10^° 


San andreas - USA 


This Paper 


1.68 


4.4 X 10^° 



movement data of the San Andreas fault (see Fisrure 1). We have also noted that although 

n 

the predicted values for q are very similar in our model and the model of Ref. |1I| , the value 
for the energy density a differ by 17 orders of magnitude, which could be used to check the 
viability of both models (a summary of main results of the present and other recent analyses 
is shown in Table I). Finally, it is worth mentioning that the estimate for the nonextensive 
parameter from the San Andreas Fault data considered in this paper is in full agreement 
with the upper limit q < 2 obtained from several independent studies involving the Tsallis 
nonextensive framework J22|. 
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